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Abstract 



We consider potential scattering theory of a nonrelativistic quantum me- 
chanical 2-particle system in R 2 with anyon statistics. Sufficient conditions 
are given which guarantee the existence of M0ller operators and the unitarity 
of the S'-matrix. As examples the rotationally invariant potential well and the 
(5-function potential are discussed in detail. In case of a general rotationally 
invariant potential the angular momentum decomposition leads to a theory of 
Jost functions. The anyon statistics parameter gives rise to an interpolation 
for angular momenta analogous to the Regge trajectories for complex angular 
momenta. Levinson's theorem is adapted to the present context. In partic- 
ular we find that in case of a zero energy resonance the statistics parameter 
can be determined from the scattering phase. 
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I. INTRODUCTION 



In recent years the theory of identical quantum mechanical particles with braid group 
statistics has received increasing attention (for a recent review and references on anyon 
physics see e.g. |35| , |39D . The original observation that in two dimensional configuration space 
IR 2 identical particles mayobey statistics different from Bose or Fermi statistics is due to J.M. 
Leinaas and J. MyrheimtJ. They based the notion of quantum statistics on the topological 
structure of the classical configuration space for identical particles. The relevant symmetry 
group is then shown to be the braid group which replaces the permutation group. Models 
for particles with a one-dimensional representation of the braid group were first discussed 
by F. Wilczek, who coined the name anyons for particles with these new statistics^! (see 
also |2H|JZ2] , |2"T| ) . Wilczek suggested the following physical picture of anyons: Magnetic flux 
tubes (vortices) are attached to either charged bosons or fermions. The latter then give rise 
to arbitrary Aharonov-Bohm phases when transported along paths exchanging the parti- 
cle positions. The magnetic flux tubes are described by long range gauge potentials whose 
curvature vanishes and can be related to Chern-Simons theory [56],[57[ This point of view 
was taken up by several authors (e.g. |7|,|23]) leading to a second quantized version of anyons 
obtained by coupling a Chern-Simons U(l)-gauge potential to a matter fieldB 
Also the general case, where the finite dimensional representation of the braid group 
is not one-dimensional, has been considered. The corresponding particles are called 
"plektons"liiJii. The relevance of braid group statistics in conformal quantum field theory 
was realized by Tsuchiya and KanieU and in algebraicquantum field theory by J. Fr6hlich@ 
and K. Fredenhagen, K.-H. Rehren and B. SchroerBH 

However, in the discussion of particles with braid group statistics the main focus is on anyons 
(abelian statistics), since in many particle theory they might provide an explanation for the 



fractional quantum Hall effect and for high T c superconductivity ( for a review see e.g. ffq.pC 
and |58|j36|) . In this article we want do discuss non-relativistic two-particle potential scat- 
tering theory for anyons. This is done in the framework of ordinary quantum mechanics, i.e. 
we insert a gauge potential of the Aharonov-Bohm type in the center of mass Hamiltonian. 
Emphasis is then put on showing how well known techniques of scattering theory extend to 
the case of anyon statistics. The motivation is twofold. On the one hand scattering theory 
is a powerful tool in spectral analysis and thus might be helpful for a better understanding 
of fractional statistics. On the other hand scattering data can be used to compute the virial 
coefficients of an interacting anyon gas!. Hence, it is possible to infer bulk properties from 
scattering theory. The latter are of main interest in the investigation of the above mentioned 
phenomena. 

The first calculation of the second virial coefficient for the two-particle anyon system was 
done by Arovas et alM They considered the case when only the statistics is present cor- 
responding to the case of Aharonov-Bohm scattering0@. See [5t] and references therein for 
recent articles on the second virial coefficient of interacting anyon systems. In a forthcoming 
publication we intend to apply the results of this article to the calculation of the second 
virial coefficient. 

The article is organized as follows. In section |IT] we review the quantum mechanics of two 
"free" anyons in order to introduce our notation and to keep the paper self-contained. 
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In particular, we give the energy eigenfunctions, the resolvent (Green's function), and the 
propagator. We also recapitulate the relation of the free anyon system and Aharonov-Bohm 
scattering which will be used in our discussion of the differential cross-section in section 
[V , In appendix [A] we recall the equivalent differential geometric formulation in terms of 



vector bundles, which in the case of anyons are line bundles. In particular, we recall that 
there exists a canonical hermitian connection encoding the statistics such that the "free" 
Hamiltonian is the canonically associated Bochner Laplacean. 

In section |TTl| we consider scattering theory for the interacting anyon system obtained by 
adding a potential to the center of mass Hamiltonian. We give sufficient conditions for the 
existence of the M0ller operators, which also cover the non spherical symmetric case. Ap- 
plying the Kuroda-Birman theoremil we derive the unitarity of the resulting S— matrix. 
In section |IV| we discuss the differential cross-section with the modifications necessary to 
accommodate anyon statistics. We show that the scattering amplitude splits into two parts, 
one describing the effects of the statistics the other the interaction represented by the po- 
tential. 

In section [VJ Jost functions are introduced which depend on the statistics parameter for 
anyons. The latter enters in the form of continuous angular momentum. This establishes a 
connection with Regge trajectories in the theory of complex angular momenta. We conclude 
section by showing how Levinson's theorem, which relates the scattering phase shift to the 
number of bound states, carries over to the present situation. In case that the Jost function 
vanishes at zero energy we derive an explicit formula giving the statistics parameter in terms 
of the scattering phase shift. 

Section [VTj and |V1J| are devoted to explicitly solvable examples. In section [VTJ we examine 
the 5-potential. This case also figures under the name of anyons without hard-core condi- 
tion. The corresponding resolvent is calculated in closed form in appendix [D] and the bound 
state problem is then considered. We also remark on the modification of Levinson's theorem 
and find an additional formula relating the statistics parameter to the scattering phase. In 
section [VII| we discuss the square well potential. The Jost function is calculated and Regge 
trajectories are plotted which show the dependence of the point spectrum on angular mo- 
mentum and the statistics parameter. In section [VTJ and |VU| we provide numerical examples 
for the differential cross-section which display the interpolation between Bose and Fermi 
statistics when the statistics parameter for anyons is varied. 

The results presented here are based in part on the diploma thesis of two of the authors 
(C.K.0 and G.lH). 

Throughout the article we will work in atomic units, h = e = m = 1, where m denotes 
the mass of the particles. In particular this sets the reduced mass of a two-particle system 
equal to 1/2. In estimates we make the convention that C, C(e) etc. denote generic constants 
depending on e etc. 



II. QUANTUM MECHANICS OF TWO ANYONS 

The theory of identical particles with statistics differing from Bose or Fermi statistics may 
show up when the configuration (or momentum) space of one particle is the two dimensional 
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Euclidean space. Henceforth we will often use the complex plane C to describe such a 
space. For the configuration space C x C of two non identical particles with points labeled 
by (-21,-22) the relative coordinate z = Z\ — z 2 changes into — z if the coordinates z\ and 
z 2 of the two particles are interchanged. The basic observation of Leinaas and Myrheim 
was, that by leaving out the case where the two particles are at the same point, i.e. where 
z — 0, the configuration space in the center of mass frame of two identical particles should 
be the space obtained from C* = C \ {0} by identifying the points z and —z. This space 
is therefore the orbit space C*/Z 2 , where Z 2 = {+1,-1} acts in the obvious way as a 
transformation group on C*. This space is also obtained from the closed upper half-plane 
HI in C minus the origin, i.e the set HI \ {0}, by identifying the points x and — x on the real 
axis. Geometrically this leads to a cone with removed apex as configuration space. The 
obvious choice of polar coordinates (r, 9) e M + x [0, n) on HI carries over to the cone C*/Z 2 . 
We take the state vectors of the system to be the square integrable functions on the cone 
£ 2 (C*/Z 2 ) or equivalently the 7r-periodic functions on the punctured plane C*. The fact 
that the cone's apex or respectively the plane's origin is removed, allows for the particles 
to carry flux-tubes. This corresponds to the physical picture introduced by F. Wilczek (see 



e.g. |57|j58| ), who named such particles anyons. The flux-tubes are taken into account by 



inserting a gauge potential of Aharonov-Bohm type 

OL 

A a = -e e , are [0,1] (II.l) 

into the center of mass Hamiltonian. Here a is the so called statistics parameter and eg 
denotes the unit vector corresponding to the polar angle. Choosing units in such a way that 
h — 1 and setting the mass of the particles to one, the resulting Hamiltonian has the form 

2 d 2 1 d 1 / d 

a V + z a ) r g r r 2 

in polar coordinates. This operator will be considered to be the free Hamiltonian in the 
center of mass frame for a two-particle anyon system with statistics parameter a e [0, 1]. If 
a = the particles actually behave like bosons, while for a = 1 they behave like fermions. 
In appendix [X] we give a short review of a mathematically precise formulation of this model 
in terms of vector bundles, following ^0|,^l] (see also ^5],^6[ There we also argue why a = 0, 1 
corresponds to bosons and fermions respectively and why we may restrict the parameter a 
to the interval [0, 1]. 

We now determine the spectrum and the eigenf unctions of H (a). We start with the 
decomposition 

£ 2 (C k /Z 2 )=£ 2 (R + ,rdr)®£ 2 (S 1 ,de) , 

where the points in S 1 , the unit circle in C, are parameterized as exp(2i9) with < 9 < ir. 
This leads to the decomposition 

£ 2 (C*/Z 2 ) = © f) 2m , l) 2m = £ 2 (M+, rdr) © <^ \ . (II.2) 

meZ V 71 



This decomposition is of course related to the following fact. The rotation group £0(2) = 
U(l) maps C* into itself and commutes with the action of Z 2 , thus acts on C*/Z 2 and defines 
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a unitary action on L 2 (C* /Z2). Also Hq(ol) commutes with this action of U(l), it is diagonal 
w.r.t. the decomposition, i.e. 

H (a) = © ^o,2m(a) 

m£Z 

with 

TT , , [ d 2 I d (2m + a) 2 
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\ o?r 2 r dr r 

on £ 2 (M + ,r<ir) = f)2m- To find all solutions of the stationary Schrddinger equation for 
H (a) it therefore suffices to find the solutions of the Bessel equation for each m G Z 



dr r dr r 



+ Z +E\R(r) = 0. (II.3) 



For definiteness we choose those solutions of equations ( |11.3|) which are regular near r = 0. 
This gives the following improper eigenfunctions 

e 2ime 

(j) a] m,E(r,9) = —= ■ J\ 2m+a \{VEr) . 

V Z7T 

In what follows we will use the notations 

k := \Z~E~, /j, :—\ 2m + a \ . 

The identity 

,00 1 

/ Jp(at)Jp(bt) tdt = —= 5 (a — b) , for a, b > 
Jo V ab 

gives the following orthogonality and completeness relations 

(4>a;m,E \ 4>a;m>,E>) = 5 m ,m' S(E - E') (II.4) 

J2 I dE <P«;™Ar, 0) K; m ,E{r',E') = -6{r- r') 5(9 - 9') . (II.5) 

We will also need the integral kernel for the resolvent i2o )Q (/Z) := (Ho( a ) — z )~~ a l so called 
Green's function, given by 

+00 

(r, (9 I Ro, a (k ± fe) | r', 9') = ±- ^ e 2 ^*-*') J M (£r<) flJ(Ar>) . (II.6) 

m=— 00 

Here we have used the standard convention r> := max(r, r') and r< := min(r, r'). To 
abbreviate notation we will sometimes make the convention that Hrr denotes the first and 

second Hankel function Hp and if^ 2 ' 1 respectively. J M is the modified Bessel function and 
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Kfj, is the MacDonald function (see e.g. [54]). To establish ( |II.6|) one uses the well-known 
formula (see e.g. |37l) . 

r kdk J ^ J ^ T,) = 7 M (cr<) K^ cr> ) [Re// > -1, Rec> 0]. 
Jo fc + c 

Finally we give the kernel of of the unitary time evolution exp (—itH (a)): 

& a {r, 6; r', 6'; t) := (r, 6 \ e - itH °^ \ r', 6') (II.7) 

We can make this explicit by adapting a calculation, which is quite standard in the context 
of the Aharonov-Bohm effect0. Along the lines of (4|,|19| we obtain the following result: 

£a(r, 9; r', 6'; t) = ^, (r, 9; r', 6'; t) + k a (r, 6; r', 6'; t) . (11.8) 
Here the two subexpressions and h a are given as 

8. a0 (r,6;r',6';t) =-L e~U^) e ta(*-*'- fag n(*-*')) x ( IL9 ) 

- — sgn(# - 0') + — cos(fl - 0'; 

k a (r, 6; r', 6'; t) =-*- ■ e^ 2 ^' 2 ) x J a (£, - 0') , (11.10) 

where we have introduced 



-ya 



I a (p, X )= dye^ - (11.11 



2y-2i x ' 

The formal relation to the Aharonov-Bohm effect used in the derivation of ( |T9 ) and (|H.10| ) 
does not come by accident. In fact, the Hamiltonian of the two-anyon system coincides with 
the Hamiltonian of the Aharonov-Bohm effect when restricted to the subspace of symmetric 
wave functions. As has been realized before (see e.g. |58|) one can exploit this by describing 
the anyonic dynamics with the help of Aharonov-Bohm scattering and thus demonstrating 
the non-trivial character of the statistical interaction. The description of the Aharonov- 
Bohm effect in terms of scattering theory was first given by Aharonov and Bohm themselves 
and later taken up by several authors, among others ^4,48,49. We will mostly follow the 
discussion presented in [45] because there the time-dependent as well as the time-independent 
scattering formalism are considered. The wave operators and the scattering operator of the 
Aharonov-Bohm effect are formally defined by 

tt± B := s -\ime itHoia) e- ltHo and S AB ■= {^abY^ab ( IL12 ) 

t— >±oo 

respectively. Here and henceforth H denotes the bosonic Hamiltonian H (a = 0) and the 
symbol s-lim stands for the strong operator limit. Furthermore, by writing ( |H.12| ) we have 
implied the restriction of the Aharonov-Bohm scattering to the subspace of symmetric func- 
tions or equivalently to the space £ 2 (C*/Z 2 ). It has been shown@ that the wave operators 
exist and are complete, whence the scattering operator is unitary. For later use we give the 
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explicit form of the integral kernels of f^ B , i.e. the stationary scattering states, which can 



be obtained by symmetrizing the results given in |48|.|49 



2m+ 



a \(kr)e 



i2m(9-0') 



(11.13) 



where \k,9') denotes the symmetric plane wave, i.e. 

(r, 9 | k, 9') = \ cos(At cos(# - 9')) . 

The scattering phase shift corresponding to Sab was first derived by Hennebergerii. Using 
his result we can rewrite the scattering operator in the compact form 

Sab = e ina P> + e - i ™P < (11.14) 

with P> and P< denoting the spectral projections onto the subspaces of positive respectively 
negative angular momentum. Equation ( |11.14j ) serves physical intuition by clarifying the 
effect of the gauge potential A a defined in (|II.1|) . 

However, of practical importance is the integral kernel of the S'-matrix in momentum space 
which was first given in |48|. After symmetrization we obtain for a e [0, 1] 



(k, 9 | S AB I k', 9') = 45(k 2 - k 



/2\ 



sin ixcx 

cos wa (5(9) + PV (1 + z cot 6) 

2ixi 



(11.15) 



Here, G = 9 — 9' and PV stands for the principal value prescription. From ( |11.15| ) one 
immediately derives the scattering amplitude 



fA B (k,e) 



(ik) 
(ik) 



[9 \T AB (E = k 2 )\9 ! ) 



(11.16) 



SIT! TTf\' 

(cos ira - 1) 5{Q) + ^— PV (1 + % cot 6) 



where {9 \ Tab(E) \ 9') denotes the on-shell matrix element of Tab '■= Sab — 1 defined by 
the equation 

(k, 9\T AB \ k', 9') = 26(E - E') (9 | T AB (E) \ 9') 

with E = k 2 and E' = k' 2 . Note, that equation ( |11.16| ) differs from what one would get when 
symmetrizing the amplitude in the way calculated by Aharonov and BohmS. In their result 
the contribution of the (^-function was left out. This violates the unitarity of the S'-matrix 
as was pointed out by RuijsenaarsS. Thus, away from the forward direction we end up with 
the following expression for the differential cross-section 

daAB \f AB {kM 2 



d9 



sin na 



iik 



(l + cot 2 #), 9y^0, 



(11.17) 



which has a nontrivial angular dependence. Note that the total cross-section is infinite if 
a ^ Z because of the singular contribution in the forward direction as displayed in ( |H.16| ). 
The latter has been interpreted to be characteristic of anyon statistics reflecting the long 
range nature of the statistical interactionil' p23 . For a = 0, 1, i.e. for bosons and fermions, 
the cross section vanishes. This is obvious from equation QH.14j) which shows that bosons 
are not scattered at all while fermions pick up a factor minus one. 
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III. POTENTIAL SCATTERING 



In this section we discuss two-particle potential scattering for anyons. Time dependent 
scattering theory involves the comparison between two dynamics, one of which is considered 
to be "free" in an appropriate sense. Here, we shall consider the "free" dynamics to be 
the one defined by H (a). This appears to be a natural choice, if one views statistics to 
be an inherent property of the particles. This is a generalization of the fermionic picture, 
where usually the free time evolution acts on the space of antisymmetric functions. Note, 
however, that this "free" dynamics contains an interaction given by the long range gauge 
forces encoding the statistics. In particular, non-integer values of a lead to a highly non- 
trivial "free" time evolution similar to the Aharonov-Bohm effect (see our discussion at the 
end of section [TT]). As a consequence not all methods used in potential scattering theory 
may be applied. For example, we do not know how two adopt Enss'El method, which makes 
use of Fourier transformation and is therefore not suitable in the present context. Thus, in 
order to prove the existence of the wave operators and the unitarity of the S'-matrix, we 
will rely on Cook's method (see e.g. 06|) and the Kuroda-Birman theoremU respectively. 



Following our discussion in section |TJ we will formulate the two particle scattering theory in 
the center of mass system using the Hilbert space C 2 (C*/Z 2 ) with H (a) = —A a being the 
free Hamiltonian. Let V, the potential, be a measurable function on C*/7*2- For example V 
may result from a function, also denoted by V(z), on C with V(—z) = V(z) and which acts 
multiplication operator. We set 

H(a) := H (a) + V = — A Q + V 

as an operator on C 2 (C* /Z2). At this point we are not concerned with giving a criterium for 
self-adjointness of Ho(a); below we specify a certain class of spherically symmetric potentials 
for which this operator is self-adjoint. The wave operators f° r the pair (H(a),Ho(a)) 
are defined as follows 

Q± := s-lim e itH(a) e - itH ° {a) 

provided the strong operator limit, denoted by s-lim, exists. Note that in the present context 
the absolute continuous spectrum of H (a) is the positive real axis including the origin and 
the associated space is all of £ 2 (C*/Z2). The S'-matrix is then defined as 

s a : = (n+)*n- 

Now we have the 

Theorem III.l. Let (l+r)V be in £ 2 (C*/Z 2 ). Then the wave operators for the pair 
(H(a), H (a)) exist for all a. 



According to Cook's theorem (see e.g. |46|) it suffices to prove the following lemma. 
Lemma III.l. Under the conditions on V stated in the Theorem one has 

\\V e~ isHo{a) (h\\ ds < 00 



for a dense set of <fi 's and for some t = t(<p) < 00. 



S 



We will provide two proofs of this lemma. The first one will use the explicit form of the 
integral kernel of e~ ttH o{ a ) given in expressions (|II.9| ), ( |H.10| ) and ( 11.11 ). The second proof 
will make use of the asymptotic form of Bessel functions near the origin and at infinity. 



1. Proof of lemma UL 1 : It suffices to consider the case a G (0,1) since the cases a = 
and a = 1 are the already well known bosonic and fermionic situation respectively. We will 
choose to be of the form r0(r, 9) = -§p^(r, 9) with $ G C^°(C*/Z 2 ). It is easy to see that 
such 0's form a dense set. 
Using (|11.8|), we write 



with 



|(e-^V)(r,0)| < F a , (r,9;t)+F a (r,9;t) 



F a ,o{r, 9- 1) 



OO / 1 7T 



o jo 



£ a , (r, 9; r, 9'; t) <j>(r', 9') r'dr'd9 



F a {r,9;t) 



OO pTV 



k a {r,9;r',9';t)(f){r',9')r'dr'd9' 



>o Jo 

We start with an estimate of F a ^. Partial integration w.r.t r' gives 



F a , (r,9;t) < - J J G Q {r,9-r',9'-t)\${r',9')\dr'd9' 

where Go satisfies an estimate of the form 

< G (r,9;r',9';t) < (1 + r) C 
uniformly in r, 9,t,a G (0, 1) and (r', 9') in the support of $. This gives 

'l + r) 



F a ,o(r,e;t) < C 



t 2 



and hence 



OO Pit 



/~i POO PIT 

\V(r,9)F afi (r,9;t)\ 2 rdrd9 <- / |(1 + r) V(r, 9)\ 2 r dr d9 
'0 Jo t Jo Jo 

uniformly in a G (0, 1). We turn to an estimate of F a . By ( 1I.10| ) 



0<F a (r,9;t) < j 

V 



OO PIT 



o Jo 



with the notation 



P=— , X = 0-9'. 
H 2t 



(III.l) 
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Adding and subtracting J a (0, x) gives 
0<F Q (r,#;t) 



< 



+ 



c 
1 
c 



OO /*7T 



\I a (p, X )-Ia(0,x)\ W,ff)\ r'dr'dO' 



JO 

OO pTV 



JO 



e-TE^+r' 2 ) / a (o, x) -^-Mr', &) dr'dO' 







dr' 



(III.2) 



(III.3) 



We need the following lemma which will be proved in Appendix |B]. 
Lemma III. 2. The quantity I a (p,x) satisfies the estimates 

\Ia(P,X)\ <C 

\I a (p, X )-I a (0,x)\ <C(e)p e 

uniformly in p and x for all < e < min(a, 1 — a). 

We use this lemma combined with a partial integration w.r.t. r' in the second term of 
the r.h.s of ( |111.2| ) to obtain 



F a (r,6;t) < C(e) 



Therefore, since r e < 1 + r we finally arrive at 



r e 1 



oo rTT 



JO 



\V(r,9)\ 2 F a (r,9;t) 



rdrdO < 



C(e) 



t 2 t 1+( 



o Jo 



(l + r) 2 \V(r,9)\ 2 rdrd9. (III.4) 



Combining (|III.1|) and ( |IIL4j ) shows that \\Ve itH ow || is integrable in t on the interval [1, oo) 
say, concluding the 1 . Proof of lemma □ 



2. Proof of lemma III. 1 : For fixed and given a G (0, 1) we use the spectral decomposition 
of H (a) to obtain an a-dependent unitary equivalence 



(III.5) 



defined by 

1p\ > tp = {lp m }meZ , ^m{E) = (<t> a ;m,E I tp) , 

such that H (a) just turns into a multiplication operator H (a): 

(H^a)ip) m (E) = (H^)ij) m (E) — E$ m (E). 
Via this isomorphism we have 

||y e -itHo(a)|| 2 = 

poo poo 

/ dE dE'e- u ( E '- E ^ m (E)v™' m '(E,E')^ m ,(E') , (III.6) 

/^^JO JO 
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where 

v™' m '{E,E') = { ( f> a . rn , E \V 2 <f> a ., m ,, E ,) . 
We now choose the following dense subspace T> C V(H (a)) in £ 



V 



G L | ^ m G C^°(IR + ) and $ m = for almost every m j . 



Correspondingly we set V = (V a ) l T>. We need the following lemma, which will be proved 
in Appendix |B[ 

Lemma III. 3. Let V satisfy the conditions of the Theorem. For arbitrary m and ml the 
functions v™ ,m (E, E') have measurable partial derivatives up to order 3 in E and E' on 
(0, oo), which are essentially bounded on compact sets. 

With this lemma at hand we now proceed as follows. We use the identity 

ue _ ( 1JM JtE 
~ [itdE 



to perform 3 partial integrations w.r.t. E in ( |III.6|) . This gives for ip E V 

\\Ve- itH °W i>\\ 2 <^Y1 j dE j dE'\dl[^ m {E)v™' m \E,E')}\ \^{E') 



which again shows that \\Ve ltH °( a > ^ | is integrable in t in the interval [l,oo) say. This 
concludes the second proof of lemma |III.1| . □ 



Now we turn to the situation, where the potential is centrally symmetric, i.e. where V 



V(r). According to 11 it is possible to define H(a) = H (a) + V as a self-adjoint operator, 
if the potential is of the form - + \ + W(r) say, where 7 and f3 are arbitrary reals, b e [0, 2) 
and W is a bounded function on R + . In general H(a) is obviously diagonal w.r.t. the 
decomposition ( [II. 2| ) such that 

H(a) = © H 2m (a) with H 2m {a) = H 0>2m (a) + V 

acting on £ 2 (IR + , rdr). This leads to a corresponding decomposition for the wave operators 
and the S"-matrix 

fi± = © fi± 2m , where Q± 2m = s-lim e ltH ^ a) e - itH °^ a) (III. 7) 



and 



£>a — © S a Om With S a 2m ~ (^a2m)*^o 2m' 



The Kuroda-Birman theoremil now leads to the following result. 
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Theorem III. 2. Let the centrally symmetric potential V satisfy 



\V(r)\ r air + 



/oo 
\V{r)\dr < oo. 



Then the wave operators fi a2m exist and are complete. In particular the S -matrices S a ,2m 
are unitary and so are the S a . 

We recall that the wave operators for an arbitrary pair (H, H ) are called complete if 
they are unitary when considered as operators from the absolutely continuous subspace of 
H to the absolutely continuous subspace of H. Note that these conditions on V in the 



central symmetric case are weaker than the conditions used in theorem [III.l 
Proof. We have to show that the operator 

\V\^{H 0>2m (a) + k 2 y l V^ 

in C 2 (R + ,rdr) has finite trace class norm H-^, where V?(r) = sgn(V(r)) ■ |V(r)| 5 . Here it 
suffices to choose any k 2 > 0. By ( [II. 6|) we have 



(#o,2m(a) + k 2 ) 1 
with fi = \2m + a | . Therefore we obtain the following a priori estimate 

i r°° 

\V\* (H 0t2m (a) + k 2 )' 1 V* < / \V(r)\ |/ M (Ar)| \K„(kr)\rdr . 

1 Jo 



(III.8) 



The following estimates for 1^, ^ can be derived from the asymptotic behavior near the 
origin and at infinity (see e.g. 0,^) 



for kr < 1 : 
for kr > 1 : 



|^(At)| < C(/x)(A;r) 



\Ukr)\ <C{p) 



\K,{kr)\ <C((x) 



-kr 



Now we choose k = 1. For this choice of k the right hand side of ( |IH.8| ) is bounded by 



CM 



\V{r)\rdr+ / \V(r)\dr 



concluding the proof of the theorem. 



□ 



IV. THE DIFFERENTIAL CROSS SECTION 

We now turn to the discussion of the two-particle scattering cross-section. The conven- 
tional approach when dealing with identical particles is to compute first the cross-section for 
distinguishable particles and in a second step to symmetrize or anti-symmetrize it in order 
to describe boson or fermion scattering. This is allowed because all observables commute 
with the projection operators onto the subspaces of symmetric and antisymmetric wave 
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functions respectively (see e.g. [51| for a detailed account on this issue). In the context of 
anyon statistics this procedure is not applicable since no corresponding projection operators 
on the anyonic Hilbert spaces are available. Hence, we will investigate the two-particle 
cross-section in the center of mass system in a way similar to the single particle case. Doing 
this one encounters an additional difficulty, namely the presence of the gauge potential A a 
encoding the statistics. The latter gives rise to a non-trivial differential cross-section even 
if no dynamical potential is present, i.e. V = 0. This should not come as a surprise since 
the "free" dynamics generated by H (a) is described by Aharonov-Bohm scattering in the 
infinite time limit (see our discussion in section [TJ). Thus, the cross-section will display both 
the statistical as well as the dynamical interaction represented by A a and V respectively. 

Let us first consider the simplest case where a = 0, i.e. the particles are bosons. Then 
the natural choice of a basis for describing the scattering is given by the symmetric plane 
waves \k,9), because the incoming asymptote is usually taken to have a sharply peaked 
momentum distribution. Moreover, we recall that the localization of some state ip at large 
times can be determined by means of its Fourier transform, 

lim [ \e- itH y\ 2 = [d9kdk\ (k,6 | V) T, (IV.l) 



where C := {(k,9) : k G M + , 9 G (a, b) C [0,7r]} denotes a cone in real space on the left 
hand side and in momentum space on the right hand side. The above identity can be found 



in several text books on scattering theory, see e.g. [44. 

If en ^ Z the flux-tube destroys translation invariance, whence the plane waves are not 
eigenstates of Ho(a). Thus, we have to look for an appropriate replacement such that the 
new basis diagonalizes H (a) and we still can form an incident wave packet with sharply 
peaked momentum distribution. It turns out that this can be achieved by making use of 
Aharonov-Bohm scattering theory. We assign to each symmetric plane wave denoted by 
\k,9) the corresponding stationary Aharonov-Bohm scattering state, that is 

\k, 9) in := Q AB \k,9) and \k, 9) out := Q+ B \k,9) (IV.2) 

The symbols stand for the Aharonov-Bohm wave operators as defined in (|II.12|) . Note 
that each of the sets { \k, 9} m } and { \ k, 0) ou *|forms a complete orthonormal system since the 
Aharonov-Bohm wave operators are unitaryca. Moreover, the basis elements are (improper) 
eigenstates of H (a) due to the intertwining relation 



Their explicit form was given in section equation (|II.13|) . Now, by construction every 



wave packet in the \k,9) out basis under the "free" anyonic time evolution exp(— itHo(a)) 
approaches the corresponding wave packet in the plane wave basis \k, 9) as t — > =foo. 

in 

Therefore, we shall refer to \k, 9)° ut as the anyonic state with incoming respectively outgoing 
(relative) momentum (k,9). Note that according to their definition the "in" states are 
transformed into the "out" states by the Aharonov-Bohm S'-matrix 



out 



(k, 9 | k', 9') in = (k, 9 | Sab I V, 9') . 
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An additional argument for the usefulness of the two bases \k,9)° ut is the following gener- 
alization of Dollard's theorem to anyonic dynamics, 

lim [ \e- itHo{a) ^\ 2 = [ d9kdk\out(k,9 \^j)\ 2 . (IV.3) 

The proof is immediate. By construction the state e~ ltH °(Vl^ B )*-ip approaches the anyonic 
state e~ ttHo<y °^ if) in norm as t — > ±00. Applying ( |1V.1| ) yields the desired relation. 
We are now prepared to consider a single scattering event of two anyons when a potential V 
is present. Denote now by ip the incoming wave function. According to ( [IV. 3| ) the probability 
that the particles are scattered into C is given by 

lim / \e- itHo{a) S a ijj\ 2 = [ d9kdk\ out (k,9 | S a ^) | 2 , (IV.4) 



where S a is the scattering operator introduced in section |TJ. On physical grounds ip is 
assumed to have a sharply peaked momentum distribution at t — —00, whence it will be 
given in the basis \k, 9) tn . Therefore, we perform the following transformation 

out (k,9 I S a i[>) = [ d9'k'dk' out (k,9 \S a \ k',9') in in (k',9' I V) 

dO'k'dk' (k, 9 I (tt+ B yS a tt AB I k', 9') in (k', (IV.5) 

where in the second step we have used the defining relation ( |1V.2| ). Thus, we are lead to 
consider the operator 

which can be identified with the S"-matrix resulting from the wave operators 

fi ± ( J ff(a), J ff ) = s-hW^ (a) e-^ . 

t— >±oo 

This can be easily verified by using the well known chain rule for wave operators (see e.g. 

H) 

= n± B . (iv.6) 

Notice that the scattering operator S 1 " 1 incorporates the statistical as well as the dynamical 
interaction while S a only takes care of the latter. This is best displayed by setting the 
potential equal to zero, 

V = =3- S a = 1, S 1 " 1 = Sab- 

Therefore, in order to accommodate the anyon statistics of the particles we shall define the 
scattering amplitude f a refering to the dynamical interaction V as follows. Denote by T 1 ® 1 
the operator 

T L ot '■= (^ab)*( s * - 1)^ab = s «' _ s ab- 
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Then the scattering amplitude is given by 



/ (M - o') ■= {lY ( d I T i°\ E = k2 ) I #) , (iv.7) 

where (6 1 | T^ l (E) \ 9') denotes the on-shell matrix element defined by the relation 

(k, 9 | | k', 9') = 25{E - E') (9 | T*°*(£) | 9') 



with E = k 2 and £" = fc /2 . From ( gVjg ) and @VJ) one can now derive the differential cross- 



section analogously to the single-particle case, e.g. |5"T||B| . Away from the forward direction 
one ends up with the expression 

^(E ,9 )=7rEf*\{e\T t a ° t (E )+T AB (E )\e )\ 2 

= \f a (ko, e-e ) + f AB (k , e - 6> )| 2 (iv.8) 

with 9 ^ 9q. Here, Eq = k%, 9q determine energy and direction of the incident particle beam 
and (9 \Tab{E) \ 9') denotes the on-shell matrix element of the operator Tab — Sab — 1 
(see section [D] equation ( 11.171 ) for the explicit expression of the Aharonov-Bohm scattering 



amplitude Jab)- Note that we have a normalization factor ir instead of 2ir in ( [IV. S ) be- 



cause the polar angle is restricted to the interval [0,7r]. As mentioned in section [IT] the total 
Aharonov-Bohm cross-section is infinite for a Z, whence the total cross-section corre- 
sponding to ( |IV.8| ) is infinite as well. However, the above differential cross-section coincides 
with the usual one for bosons or fermions if we set a = and a = 1 respectively. This 
can be most easily seen by use of the defining relation ( [IV. 2|) for the "in" and "out" states 
which for the values a = 0, 1 become symmetric and antisymmetric plane waves respectively 
(compare ( |11.13| )). In case a = 1, however, this is only true up to an angular dependent 
phase factor which comes in by representing fermions as bosons with attached flux-tubes. 
This does not influence the outcome since the differential cross-section is given by the square 
modulus of the scattering amplitude. Thus, as special cases we obtain the familiar result 
that the differential cross-section for bosons and fermions is given by the square modulus 
of the symmetrized and anti-symmetrized scattering amplitude respectively. In particular, 
Jab vanishes and the total cross-section becomes finite for suitable short range interactions 
V. 



V. JOST FUNCTIONS AND LEVINSON'S THEOREM 

In this section we will introduce Jost functions^ indexed by a continuous angular mo- 
mentum and discuss their properties. This continuous parameter leads to an alternative 
formulation of the generalized Levinson theorem. We start by adapting the standard theory 
and results of Jost functions to the present situation (see e.g. ^2| , |i3| ,|3| j5T| , |46| ) . Conditions on 
the spherically symmetric potential will be presented at the appropriate places. 
Consider the unitary map U : £ 2 (R + , rdr) — > £ 2 (R + , dr) given by 

ip(r) i — > yfripir) . (V.l) 
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Under this map the Hamiltonian if 2m (a) = if ,2m(a) + V turns into the operator 

K = h ^ + V ( V -2) 

with fi = | 2m + a \ > 0, where 



In this section /i will be allowed to be any positive number. Also in this section the notation 
4>o(r; k, /i) will be reserved for the regular solutions (so called because of their behavior near 
r = 0) of the free Schrodinger equation (/i 0jAt — k 2 ) tp = given as 



Mr;k,v) = y/^J^kr) . (V.4) 

With the help of these solutions the orthogonality and completeness relations ( |II.4| ) and 
( |11.5| ) are now reformulated as 

4>o(r; k, n) <f) (r ; k,fi)dk = - 8(r - r) (V.5a) 

0o(^; fe, /x) 4> {r; k , n)dr = - 5(k - A/) . (V.5b) 
o 1 

We will also need the irregular, free solutions Xo given as 



o 



With help of these solutions the free Green's functions read (compare ( |11.6| )) 

Gq(t, r'; k, fi) = (r | (/i 0l/L4 - A; 2 =F i^Y 1 \ r') = \ M r <~, k, £t) Xo ( r >; fc > A*) • ( v - 6 ) 

Definition V.l. For given V the function <fi = <f)(r;k,fi) is the regular solution of the 
equation {h^ — k 2 ) ip = which for r — > approximates the free, regular solution 4> Q : 



As in the standard theory (see e.g. |42| , ^3|) one establishes the following facts. 
(AO) is the solution to the integral equation 

if, = 0o - G<Vif> with G<(r,r';k,fi) = Q(r-r')g (r,r';k,fi) . (V.7) 
Here 9 denotes the Heaviside step function and g is given as 

9o(r, r'; k, fi) = |[0 o (r; k, fi) Xo( r '> k, //) - x (r; A;, yu) o (r'; k,fi) \ . (V.8) 
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(BO) If the potential V satisfies the condition 

rdr\V(r)\; < oo 



o 



then the regular solution to the eigenvalue problem (h 0tU + XV — k 2 ) if) = has a 
convergent power series expansion in A of the form 

oo 

= ^ A n n with (fi n=0 = O and (fi n = -G$V<f> n -i, n > 1 . (V.9) 

n=0 

(CO) For fixed // > 0, r > and real V the solution has an analytic continuation in k 
into the complex plane with a cut along the negative imaginary axis. There one has 
the relation ( k > 0) 

0(r; - 0, /i) = e in ^ + ^ 0(r; ik, /i) = e 2i7r ^+D 0(r; -ifc + 0, //) . (V.10) 

The proofs of (AO), (B0) and (CO) are as in ordinary 3-dimensional Schrodinger theory (see 
e.g. |43|). As a byproduct of the proof one also has the two estimates: 





k 


r 


1 + 


k 


r 



'* C(fi) n+1 



dr ,r'\V{r')\ 



(V.lla) 



1 + \k\ r' 

(r;*,/i)l <CV)e |Imfc|r f-^-V +2 (V.llb) 



We will prove these estimates in the Appendix |C|. Relation (|V.10|) is a consequence of the 
relation 

We now introduce the irregular solutions of the eigenvalue equation with potential as those 
solutions which approximate the free irregular solutions Xo ^ r = °°- 

Definition V. 2. The functions x^ = X ± { r \ k, fi) are the solutions of the eigenvalue problem 
( h a — k 2 ) if) = satisfying the boundary conditions at infinity of the form 

lim e Ti{kr ~^ fl+ ^x ± (r;k,fi) = 1 . 

) — »oo 

These solutions will be called Jost solutions. With the help of these solutions and x ± 
the Green's function for the full Hamiltonian now takes a form analogous to the free case 
(see O) 

G (r,r;k,fi) = (r\(h a -k =F ze) \r) = w ± ^ , (V.12) 

where W is the Wronskian, i.e. 

Wiip!,^)^) = ifii{r) d r if) 2 (r) - ifi 2 (r) 9 r -0i(r) . 

We recall that since and x ± are solutions of the same eigenvalue equation, the Wronskian 
in ( |V.12|) is actually independent of r, whence the notation. 



Again one may establish the following facts (see e.g. |43| ) 
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(A±) The functions x ± are the solutions to the integral equations 

ip = Xo ~ G o Vi P » G o( r i r> i '■= @ ( r ' ~ r )9o(r,r'; k,fi) . 
Here, go denotes the function defined in (|V.8| ) which can be rewritten as 

9o(r, r'; k, n) = ^ [Xo ( r 5 A t )Xo ( r '; fc > A*) ~ ( r ! fc > aOXo ( r '; fc > A 4 )] 
(B±) If the potential V satisfies the condition 

POO 

/ rdr (1 + r)\V(r)\ < oo , 
Jo 



(V.13) 



(V.14) 



the irregular solutions x ± to the eigenvalue problem (ho,^ + AV^ — k 2 ) ip — with 
± Im fc > 0, ^ have a convergent power series expansion in A of the form 



X 



± 



With X^= = Xo aIld Xn = - G V Xn-l, U > 1 ■ 



n=0 



(C±) For V real and \x > 0, r > the functions x 1 * 1 are analytic in k for ± Im k > with 
continuous extensions to the real axis except possibly for a singularity of order \i — | 
at = 0. In case V satisfies 



rdr\V(r)\e vr < oo, r] > 0, 



(V.15) 



then the functions x have an analytic continuation into the domains given by 
±Im/c > — f f° r all > , r > with the exception of a branch cut from k = to 
= =f«tJ . There one has the relation ( k > 0) 

X ± (t; =f(z/c + 0),/i) = x ± ( r ; =F(^ - 0),/i) + 2z cos7r/i x ± ( r ; A 4 ) • 

The proof of these statements is similar to the one in the regular case and involves a proof 
of the bounds 



b£(r; *,/*)! <e Tlm(fe)r 



\k\ r 



x 



X 

n\ 



and 



|x ± (r;fc, y u)| < C^e™^ 



^ I |Imfc| r'=Flm(A;)r r/ l^( r I 

1 + \k\r' 





k 


r 


1 + 


k 


r 



(V.16) 



(V.17) 



for all fi > and r > 0. We will establish these bounds in appendix 0. Due to the uniqueness 
of the solutions x ± an d the multivaluedness of the Hankel functions (see e.g. |37|) one has 
the relation 



X ± (r;k,fi) = x T (r;k,n)= =p e ±l7T ^ x ± (r, -k, (j) 
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We are now in the position to introduce the Jost functions in the present context and 
establish their properties. Since x + an d X~ form a basis of solutions, is a linear combination 
of these two solutions. As in the standard theory (see e.g. £|3]) it turns out that this linear 
combination is of the form 

0(r; k, (i) = i [F(k, n)x~{r; k, fi) - F(-k, fi)x + {r; k, n)] , 

where the coefficient function F = F(k, fx) is called the Jost function given by the Wronskian 

F = \W{ X + A) • (V.18) 

The Jost function was first introduced by and named after R. Jost. Note, however, that our 
definition of F follows the convention introduced by R.G. Newton@, which differs from the 
original one given in E9|. We have the first main result concerning these Jost functions. 



Theorem V.l. Let the potential V satisfy the condition (|V.14|) . Then for fixed fi > 
F(k, jj) extends to an analytic function in k in the open upper half plane, which is continuous 
up to the real axis with the exception of the origin. 

In case the potential satisfies the stronger condition ( |V.15| ) ; F(k,fi) is analytic in k in 
{k | Imk > — |} except for a cut on the negative imaginary axis. Then F also satisfies 
the relation (k > 0) 

F(-ik-0,fj) = -e 2i7r ^F(-ik + 0,fi) + 2 cos tt/j, F(ik, /i) . 

The proof follows from the representation for F 

F{k,n) 



(V.19) 



the estimate flV.llbp for and the estimate (|V.17|) which is also valid for xt ■ To establish 
this relation we use the integral equation ( |V.7| ) for which for large r gives 



<f>{r;k,ii) 



r— »oo 1 

" 2 



X (r; k,fi)\l + - (xo (• ]k,fi)\ V<j>{- ;k,/i)) 



1 



Xo(r; k,n)\l + ~ (xo (• ; k , I*) I V(f>(- ; k, //)> 



. (V.20) 



Since x ± ( r ! approaches Xo{ r 'i for r — > oo, relation (|V.19| ) follows from this behav- 
ior. 



Note that making use of the analytic properties just proven there is an equivalent definition 
of the Jost function in terms of Fredholm theory, see e.g. Consider the operator 

Gq (k,fi)V with Gq the free Green's operator and V satisfying relation (|V14| ). Then for k 
positive imaginary (k G zK + ) the operator G^k, jj)V is trace class, whence its Fredholm 
determinant det(l + G^{k, fJ>)V) is well defined. It turns out that the Jost function is just 
the analytic continuation of this determinant to the upper half complex plane. 
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We note that the series expansion ( |V.9| ) for the regular solution gives the series expansion 

- oo 

F = l + -Y^F n with F n (k,fi) = ( Xo (-;k,n)\Vr(-;k,v)) ■ 



k 
o 



The bound ( |V.lla| ) for <p n combined with the bound (|V.17|) for Xo gives the bound 

\k\ r 



/~in poo 

\F n {k,fi)\ <C{fi)^— / f W |Imfc| - Imfe)r - 
n - Jo 1 

while the bound flV.llbD gives 



\k\ r 



\V(r)\ 



F(k } fi)-l\< C{n) / dr e d ImA; l- ImA; > - 

Jo 1 



r\V{r)\ 
+ \k\r ' 



The importance of the Jost function stems from its role in the discussion of the S'-matrix. 
Let fl^ be the wave operators for the pair (h^, h ^) such that one has (see (pi.7|) and ( |V.1| )) 

fij = un^u- 1 and s, = (n+y n- = us^u- 1 . 

Since commutes with h Q ^ there is a decomposition of in the form 



= dk S(k, fj,) 
J® 



with respect to the spectral decomposition of h 0ifl . Since the spectrum of h ^ is not degen- 
erate, S(k, fi) acts in a one-dimensional Hilbert space and therefore is a complex number of 
absolute value one, i.e. 

S(k, /i) = e m ^ , k > , /I > 
with 5(k,n) being the phase shift. We define by 

2i<5(fc,/t) _ i 

f(k,fi) : = 



\niik 

the partial wave amplitude for the angular momentum ±/i, such that we obtain the following 
partial wave decomposition of the scattering amplitude given in ( [IV. 7|) 

/«(m- 00 = \T t :\e)\e>) 

— (4L)3 I e -iTr\2m+a\ ^ e 2iS(k,\2m+a\) _ -Q ^i,2m(6-6') 

= e -i 7 r|2m+a|j( A . j | 2m + e i2m(0-0'). ( V .21) 

Here we have used the relation (k, 6 | T*°* | k', 0') = out (k, 9 \S a — l\ k', 6') in as well as the 
identities 
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\<ra;m,E 

Now set 



S^c,^) = 5 mm ,5{E - E') e 2 ^ 2m+a \\ 



By construction one has a solution of the eigenvalue equation {h tM — k 2 )ip = and (p ± satisfies 
the Lippmann-Schwinger equation 

± (- ;fc,/x) = <M" ;A;,/i)-(^-A; 2 ±2O)- 1 \/0 ± (- ; fc, //) . 

With the help of ± one may express the partial wave amplitude as 

f(k,n) = -^(f) 1 {M-;k^)\Vc) ) + (--k^)) . 

The solution + of the Lippmann-Schwinger equation has the following asymptotic behavior 
for r — > oo. 

<j>+{r;k,fx) > Mr;k, //) + \V^k f(k, //) e^-f") , 



2 

where we have used Xo( r jk,fx) > e l ^ kr ~2^ + z). Alternatively, using the relation </> 

r— >oo 

t ( Xo - Xo ) one § ets 



With the same arguments as in the standard theory (see e.g. H2,H3H3) one proves 



Lemma V.l. Let fi > be fixed and assume that V satisfies the condition ( V.14|) . Then 
one has 

(i) For k on the real axis (k G M + ) the following identities are valid: 

(f>(r]k,fi) = F(k,fi) (j) + (r;k,fi) (V.22) 

c (k ..) _ e 2i&(k^) _ F {k,fi) _ F{-k,n) , . 

H " 6 ~F(k,») ~ F{k,n) (V - 23) 

Furthermore, F(k, //) can on/?/ vanish at the origin. 

(ii) The zeros of F(k,fi) in the upper half plane Irak > all lie on the positive imaginary 
axis and are simple. In particular k G zM + is a zero of F(k, /i) iff k 2 is the energy of 
a bound state of with angular momentum ±/x. 

(Hi) The following limit relation is valid in the closed upper half plane 

lim F(k,fi) = 1 . 
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Since the proof is analogous to the one given in three dimensional Schrodinger 
theoryHHUll we omit it here. Now, the famous Levinson's theorem follows as a corollary by 
use of the residue theorem. Only the so called resonance case, i.e. F(k = 0,/i) = 0, requires 
some more work. 



Theorem V.2. Let /i > and V satisfy ( |V.14| ). If F(k = 0,/i) 7^ one has the following 



relation between the phase shift 5(k, /i) and the number of bound states with angular 
momentum ±/i ; 

6(k = 0, /i) - 5(k = 00, n) = 7m M . (V.24) 

Proof. The proof of equation ( |V.24j ) is standard and can be found in e.g. |42",|43,46l,51. For 
sake of completeness we recall the main arguments. Integrate the logarithmic derivative of 
F with respect to k over a closed semi-circle in the complex upper half plane centered at 
the origin. According to the preceding lemma the integrand has simple poles on the positive 
imaginary axis each of which corresponds to a bound state of h^. While the contribution 
of the semi-circle vanishes in the limit of infinite radius (see ( |U1.4| )) the integration over 
the real axis yields the phase difference by equation ( |V.23|) . Applying the residue theorem 



completes the proof. □ 

We comment on the results of the above theorem. Equation ( |III.5|) generalizes Levinson's 
theorem to continuous angular momentum /1. One might wonder what happens when fi 
varies. On physical grounds one expects that for higher angular momentum /1 there are 
fewer bound states since then the centrifugal barrier is stronger. Thus, when /i is steadily 
increased the zeros of the Jost function F move along the positive imaginary axis towards 
the origin. Every time a zero moves out of the upper half plane we get a discontinuous jump 
in the phase shift indicating the disappearance of a bound state. 

Theorem V.3. Let (jL Z, /i > and assume V satisfies (|V.15|) . In case F(k = 0,/i) = 
the above form of Levinson's theorem remains valid when fi > 1 and is modified as 

8{k = 0,/i) - 5(k = oo,//) = irfan + fi) (V.25) 

when /i < 1. 

So far we have not been able to settle the case /16Z corresponding to boson and fermion 
statistics. For real anyons, however, note that in case of positive angular momentum we have 
/i = a when /i < 1. Thus, equation (|V.25|) then specializes to 



S(k = 0, a) — S(k = 00, a) = n^n^ + a), (V.26) 

showing that one can determine the statistics parameter from the scattering phase inde- 
pendently of the detailed form of the short range potential V. Below we will find a similar 
formula for the delta potential. Setting a = 1/2 (i.e. for semions) in (|V.25|) we obtain the 
well known resonance case of three dimensional Schr0dinger theory^, since then the radial 
Schr0dinger equation is formally identical to the one in three dimensions (compare ( |V.2|) 
and ()V3D). 
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Proof. The strategy for proving (|V.25 ) is similar to the one used for ( |V.24| ) . The techniques 
used are analogous to those in three dimensional Schrodinger theory (see £|^,[5l|) but the line 
of argument is slightly changed. 

Because F(k = 0,/x) = we change the integration contour as follows. We cut the semi- 
circle at the origin and insert a small semi-circle C p of radius p <C 1 centered at the origin. 
Then there is an additional contribution to equation (|V.24j) namely, 



oo,//) + lim — 
p-+o zi 



(V.27) 



Co 



where n'^ is the number of bound states with the exception of a possible zero energy eigen- 
state. In order to compute the second term of the right hand side in the above equation 
we only need to know the small energy behavior of F. Since we have assumed that the 
potential is decaying exponentially we can continue F in the lower half plane and expand it 
into a convergent power series around the origin for suitably small k. We start by rewriting 
expression (|V.19|) , 



F(k, n) = 1 + \ (x (• 5 k, (i.) I V<j)(- ;k,n)) 



1 + 



i 



k sin 7r^t 



o(- ! k, -n) | V(j>{- ; k,/i)) 

/ ]k,fi) | V(j)(- ;k,fi 



-tTTfJ, 



(V.28) 



Here we have used the identity H, 



i 



I sm 7T/t 



[J- 



e t7Tfl Jfi{z)\ which can be found 



in e.g. |5^ , |37] . Note that in case of integer /i the Hankel function is defined by the limit of 
r.h.s. of the last equation. Then the above decomposition (|V.28| ) is no longer valid, whence 
we restrict ourselves to the non-integer case. Both integrals in ( |V.28| ) converge uniformly as 
can be shown by use of the estimates 



\(f>o(r;k,n)(f>{r;k,ii)\ < C{p)e 



2| Imfc| r 
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\Mr;k,-fi)<P(r;k,p)\ <C(p)e 2 ^ r 
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which follow from ( |V.lla|) and ( |V.llb| ). The functions (fio(r; k, —(/,) and 4>o(r; k,fi) can be 



written as a convergent power series in k 2 times a factor k~^ + ^ and k^ + ^ respectively. 
This follows from their definition ( |V.4| ) and the series expansions of the Bessel functions 
J-v(kr), J^(kr), see e.g. |54l , |3~7] From ([V.9D we see that also 0(r; k,fi) can be written as a 
convergent power series in k 2 times /c^a, since the Green's function is an even function in 
k. Therefore, in a suitably small chosen neighborhood of k = the Jost function can be 
expanded in a convergent series of the form 

oo oo 

F(k, fi) = J2 a n{v)k 2n + k 2 ^ J2 b n {p)k 2n , \k\ < 1. 



(V.29) 



n=0 



n=0 



Recall our convention to place te branch cut of F along the negative imaginary axis. The 
assumption F(0,fi) = gives ao = 0. In order to deduce from this expansion the small 
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energy behavior of F we need to know the first non- vanishing coefficient in ( V.29|) . We 
claim ai(fjj) for /x > 1 and bo (/x) 7^ for /x < 1. But then we have 

^ ^ - / °(^) if ^ < 1 
- \ Oik 2 ) if /x> 1 ' 

Thus, evaluating the integral over the semi-circle C p yields 

hrnl/ d\nF=\-^ . (V.30) 

p-*0 2iJ Cp \ -7T if /X > 1 

In order to prove the claim we consider the limit lim^ k~ 1 ~§kF. If we can show that 

then the claim follows as can be infered directly from ( V.29| ). We calculate the expression 
JgF by making use of the identity ( |V.18| ). Since the solutions x + , <\> are ill defined for k — ► 
we consider instead the modified solutions 

^(r; k, fi) := k^~^x + i r 'i k, fi) and $(r; A;, /x) := £T M ~20( r ; £;, /1) 

which are both finite and non-vanishing at k = 0. This can be seen by taking the limit 
k — > in the defining integral equations (|V.7|) and ( |V.13| ) respectively. In particular one 
then derives that <3>(r; 0, /x) behaves like for r -> and \I/(r; 0, /x) like r~ M+ 2 for r — > 00. 
In terms of $, \l/ the Jost function and its derivative w.r.t. k (denoted by a dot) then read 

F = W(ty,$) and F = + 

Since $ is a solution of the radial equation one easily verifies the identity 
d 

— W[V{r; k, /x), #(r; fc', /x)] = (k 2 - k' 2 )^(r; k, /x)#(r; fc', /x). 
Differentiating w.r.t. and setting subsequently k = k! one obtains from this the relation 

/>oo 

W[¥(r; fc, /x), #(r; k, //)] = -2Jfe / dr'#(r'; jfe, /i) 2 , (V.31) 

where the r.h.s. is finite for Imk > 0. The assumption F(0, /x) = implies that the regular 
and the irregular solution are proportional when k — 0, i.e. \l/(r; 0, /x) = ft(/i)$(r; 0, /x) holds. 
Note that the proportionality factor k(/x) is nonzero and finite because both solutions are 
non-vanishing. Thus, in the limit k — ► we have 

F(0, n) = k{^)- 1 W (*(r; 0, /x), *(r; 0, //)) + k(/x) W ($(r; 0, /x), $(r; 0, /x)). 

However, F(0, /x) does not depend on r whence we can set r = in the last equation. Since 
the regular solution $ vanishes at r = we obtain from ( |V.31| ), 
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d 

lim k- 1 —F(k, y) = lim fc -1 /^) W [tf (0; fc, /i), ^(0; k, //)] 

fc^O C/K fc— >0 



oc 



-2k(^) -1 / rfr^(r;0,/i) 2 
Jo 

POO 

-2n(ji) / dr$(r;0,/i) 2 (V.32) 



o 



The irregular solution \I/ falls off like r M+ a for large r and is proportional to the regular 
solution $ at k = 0, therefore the r.h.s. of ( |V.32| ) is finite for /i > 1. This means that the 



regular solution is square integrable and hence a proper eigensolution of the full Hamiltonian 
h^. The number of bound states is thus = n' + 1, whence we read of from ( |V.27 ) and 
( |V.3CI| ) that Levinson's theorem also holds in case F(0,fi) = when fi > 1. 



For \i < 1, however, the regular solution is no longer square integrable and the expression 
( |V.32| ) is divergent. Hence the number of bound states is given by n M = n'^ and from (|V.27|) 



and (|V.30D we infer the modified equation (|V.25|) . This completes the proof. □ 



VI. THE 5 POTENTIAL 



In this section we will introduce and discuss the ^-potential, often also called the con- 
tact potential or zero-range interaction potential because it describes an interaction which 
is "non-vanishing" only if the two particles are at the same place, i.e. if z re \ = 0. In 
the literature this figures also under the notion of anyons without hard-core condition^. 
See e.g. |38| , p!0| j5| for previous articles on this subject also in connection with field theoretic 



considerations. 

To introduce the (^-potential we follow a standard strategy (see e.g. 0), i.e. we consider 
the free Hamiltonian (in the center of mass system) restricted to a definition domain of 
wave functions which vanish at z re i = 0. The resulting operator has deficiency indices (1,1) 
and so there is a one-parameter family of self-adjoint extensions, one of which is of course 
the ordinary free Hamiltonian given as the Friedrich's extension. We discuss the bound 
state problem and the resulting scattering theory. We also check the validity of Levinson's 
theorem. 

For a start, consider the operator 

^2 _ 4 -l 

h ,ti = -d r H — 2 

on £ 2 (W + ,dr) (recall ( |V.lp ) with domain of definition T>(h ^) consisting of wave functions 
ip{r) having compact support away from the origin such that ip and d r ip are locally absolutely 
continuous and such that h ^ip (defined in the sense of distributions) is an element of 
£ 2 (]R + , dr). Recall that locally absolutely continuous functions are such that their derivatives 
in the sense of distributions are locally integrable functions (see e.g. 52). 



Let us now investigate the most general solution of the equation ho nip = 0, given by 



i 



C\ ■ r2 T ' 1 + c 2 • r 2 M 
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Using Weyl's terminology (consult e.g. [55]), one has the limit point case at r = oo. At r = 
one also has the limit point case if and only if /x > 1, otherwise one has the limit circle case. 
If now < ji < 1, by Weyl's alternative there is exactly one solution of h ^ip — 0, which is 
C? near infinity. By construction this solution must be C 2 at both infinity and at zero, and 
thus is follows that dim(ker(/io lAt ± i)) = 1 by Weyl's alternative. Therefore the deficiency 
indices are (1,1) and there is a one-parameter family of s.a. extensions (see e.g. |47|) . 

In terms of the statistics parameter a we have the following situation: If one chooses 
a = 1, then fi = \2m + a\ > 1 for all m G Z and all the {h ^} are essentially self- 
adjoint. Now, if a G [0, 1) we have \i < 1 if and only if m — 0. In this case, just as in the 
three dimensional case, the (^-potential affects only one of the angular momentum channels, 
namely the one associated with the smallest eigenvalue in the sense of the modulus (in three 
dimensions called the s-channel). 

Let a be in [0, 1). To construct the s.a. extensions of the operators h a with domain 
V(h a ) given above, following we define the regular and irregular solutions of h 0ja ip = 0: 

F a (r) := r2 +a (regular) 

pro 

G a (r) := F a (r) dr' {F a {r'))- 2 (irregular) 

J r 

where r G M + may be chosen arbitrary. The choice r = 1 gives 



G a {r) 



— -y/rlnr a = 

In order to explicitly construct the s.a. extensions of the {h 0;a } we take recourse to the 



following theorem, to be found in e.g. 55 



Theorem VI. 1. For a < 1 the operator ho yCe with domain T>{ho^ a ) has a one-parameter 
family of s.a. extensions {/io,a(<s)} s eR. which are essentially s.a. on the domains 

% a W) = W^W I limW(Va )S ,<0)(r) = O}. 

rj.0 

Here again W is the Wronskian and 

The choice s = oo gives the Friedrich's extension. 
In appendix we will prove the following 

Theorem VI. 2. The integral kernel of the resolvent of h a (s) for a < 1 is given as 
1 



ho,a( s ) - k 2 



[r, r') 



yVrr 7 J a {k r< )H^{k r> ) - A{k,a;s) ■ \fr~r~' H^\kr) H^\kr') , (VI.2) 
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where k satisfies < axg(k) < n. A(k,a; s) is given as: 



71 



A(k,a;s) = - 



-2a 



(2sa 



1 r(l + a) 

1 ) — r + COt TXCi 



sin Tia r(l 



a] 



(VI.3) 



Note that lim^oo A(k, a; s) = such that with ( |V1.2j ) in the limit s — > oo we indeed 
recover the kernel of the resolvent of the Friedrich's extension ( |V.6| ). Taking the "bosonic 
limit" a — > one arrives at the expression: 



7T 



A(k,0;s) = — 



In — + 7 + s 
Li 



(7: Euler's constant) 



Combined with ([VI. 2|) we have arrived at an expression for the resolvent of the ^-potential 
in the bosonic case, which agrees with ||. 

We turn to the discussion of the bound states. It can be shown that the essential spectrum 
as a set is equal to IRq and that R + is purely absolutely continuous. We therefore have to 
look for poles of the resolvent at energies Eb = k\ < 0. By (|VI.2|) we have to search for the 
poles of A(k, a; s) as a function of k and by (|VI.3| ) this gives the condition 



k 



2r/, _ 2a i-ira mod + a ) 

6 T(l -a) 



2as). 



(VIA) 



The restriction < arg(fc) < tt now implies < arg(/c 2ct ) < 27ta. Applying this condition to 
( |V1.4| ) leads to three different cases: 



(i) 1 — 2as > 0: Then ( [VI. 4|) may be solved to give 

v i-f \ .rr(i + a ) 

kb = kb(a, s) = 2i 



ni 



a 



[1 - 2as) 



1 

2o 



i.e. kb is purely imaginary. The case kb = occurs when a | 1. Note that this requires 
< s < Hence there is exactly one bound state. 

(ii) 1 — 2as < 0: Since by definition < a < 1, it is easy to see that in this case there is 
no solution to equation ( |VI.4D , only in the limit a f 1 one obtains kb = 0. This means 
that there is no bound state. 

(iii) 1 — 2as = 0: This immediately gives kb = 0, which implies that there is no bound 
state. 

For comparison we recall that in the three dimensional case the attractive ^-potential 
supports at most one bound state. A bound state exists if the particles obey bosonic 
statistics. For fermions there is none. In the present context a bound state always exists 
for a = (bosons) and disappears at the latest when a = 1 (fermions). The fact that for 
a = there is always a bound state suggests that in two dimensions there is no such thing 
as a repulsive 5-potential. 



27 



Below we argue, that case (iii) corresponds to a zero energy resonance in three dimen- 
sional Schrodinger theory. Figure [l] gives a plot of \kb(a, s)\ as a function of s for various 
values of a in case (i), i.e. when 1 — 2as > 0. To calculate the bound-state wave function 
ipb when 1 — 2as < 0, we note that 

1 1 



h , a (s)-k 2 k 2 -k 2 

for k ~ k b where P b is the orthogonal projector onto the eigenspace of energy E b = kf < 0. 
Since this pole of the resolvent is isolated, P b may be calculated as a residue. The calculation 
is easy and shows that P b is indeed one-dimensional and that the associated (normalized) 
eigenfunction ip b is given as 

i 

, , . / 2 sin net \ 2 I . . 

?p b (r;a)=[ y/r \E b \ K a (\k b \r) 

\ na J 

for a < 1 and hence \k b (a,s)\ > 0. For any fixed a < 1 the wave function decays expo- 
nentially for large r, which guarantees square-integrability. In the limit a | 1 the bound 
state wave function approaches r~z times a normalization constant, which is not square 
integrable. The case a = 1 should therefore be considered as a resonance at E = 0. 

We turn to a discussion of the resulting scattering theory. First we note that due to 
relation (|V1.2| ) the difference of the resolvents for ho, a (s) and the free Hamiltonian ho, a {s = 
oo) is a rank-one operator and therefore in particular trace class. Therefore the S'-matrix 
exists and is unitary by the Kuroda-Birman theorem (see e.g. ^). 

Theorem VI. 3. For given a and E = k 2 > the outgoing states for the Hamiltonian 
ho, a (s) ar e given as 

2i 

(f) + (r; k, a) = (p (r; k, a) + —A(k, a; s) ■ Xo{r; k, a) . 

7T 

In particular the resulting partial wave amplitude and the phase shift take the form 

f{k } a; s) = r=r - , (VI.5a) 

Aj 

e 2i6(k,a;s) = j + _A(k, d' s) . (VI.5b) 
7T 

Proof. We start by rewriting the resolvent ( |V1.2| ) as 

1 / 2i \ i 

- — _ _( r , r') = </>o(r<; k, a) + —A(k, a; s)xo (r<; k, a) T xt( r >i k > a ) ■ ( VL6 ) 

ho, a {s) - k 2 \ ix J k 

Now we observe that (see e.g. 0, p. 37) 

lim lim e -^+iey ( r y\ (vi. 7) 

elOr'^oo h^ a (s) - [k + IE) 2 

must be proportional to the outgoing solution + (r; k, a). The claim now follows by applying 
( |VL7D to ED . □ 
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In view of Levinson's theorem we would like to make the following comments on this last 
result: The limiting values of 8(k,a;s) at k — > + and k — > oo are a priori determined 
up to an additive constant in tt ■ Z, only. We shall choose these constants in such a way, 
that S(k, a; s) as a function of k G M + can be defined as a continuous function of k. Now, 
exp(i2<5(/c, a; s)) is of the form (a(k, a; s)+i)(a(k, a; s) — Here a(fc, a; s) is a real-valued, 
continuous function of k. Explicitly 

a(k, a; s) = ~^A(k, a; s)^ 1 + i 

(This observation also shows that the right hand side of ( |V1.5bj ) is indeed unimodular.) 
Hence 

8(k, a; s) = axg(a(k, a; s) + i) . 

Here the branch of arg(-) is chosen such that S(k, a; s) is a continuous function of all its 
arguments. In particular we can achieve S(k, a; s) G (0, n), for any k and fixed a < 1 and s. 
For the cases, where there is a bound state, this is consistent with the attractive nature of 
the interaction. By the explicit form of a(k, a; s) it is easily seen, that 

jf/, x fc^o+ 1 - sgn(2as - 1) 
d(k,a;s) ► vr 

5(k,a;s) fc ^°°> rra 

The fraction n := 1 ~ sgn ^ as ~ 1 ) takes values in {0, 1} and is in fact nothing but the number 
of bound states given by ( |V1.4| ). We therefore have the following relation 

<5(0 + , a; s) — 8(oo, a; s) = 7r(n — a) , (VI. 8) 

which is not of the form one would naively expect in view Levinson's theorem |V.2| . In this 
context we note, that the statement of Levinson's theorem also fails in the non-anyonic three 
dimensional case, if a ^-interaction is considered. However, we obtain a formula similar to 
( |V.26| ) relating the phase shift to the statistics parameter a. 



By relations ( [IV. 8|) and ( V.21[) , we obtain the following expression for the differential cross- 
section: 

^= \f(k,a;s)-e- i ™ + f AB (k,e)\ 2 , 



where Jab and / are given by expression ( II.16Q and ( VI.5a|) , respectively. Inserting these 



explicit expressions, one obtains — away from the forward direction — the following relation 



da 1 
d9 irk 



sin Tia (cot 6 - i) + e~ l7TOl —A(k, a; s) 

7T 



2 



So far, we have been unable to analyze this in a non-numerical way. Figures and |3| show 
some numerical results of the differential cross-section as a function of the scattering angle 
9. Displayed is the deviation from pure Aharonov-Bohm scattering which is asymptotically 



29 



reached when s — > ±00. Note the simultaneous cross over at 9 = na for varying s in each 
of the plots with s 7^ ±00. In particular, the point corresponding to the cross over always 
lies on the graph (s = ±00) associated with the Aharonov-Bohm cross-section. We like 
to point out that in both figures || and [| only for s = 0, 2 and a = 0.5 (semions) the 
differential cross-section actually vanishes for some scattering angle 9. In the other cases 
the cross-section becomes very small but is non- vanishing for all 9. 



VII. THE SQUARE WELL POTENTIAL 

Another explicitly solvable problem in scattering theory is given by the square well, 

"M = {o?'r>2 With V °>°- 

The regular solution of the eigenvalue problem (h^ — k 2 )ip = in the interval < r < d is 
given by 

4>(r;k,n) = 2 4> (r;q,ij), where q= y/k 2 + V . 

For r > d the eigensolutions of obey the free radial equation and hence the Jost solutions 
in this interval are just given by the free solutions, 

X ± (r;k,fi) = Xo (r;k,n), r > d. 

Making the usual requirement that all solutions are continuously differentiable we can eval- 
uate the Wronskian W( y x ± ^ ( P) at r = d and obtain the following expression for the Jost 
function, 

F ( k ^) = (|) 2 Xo(d]k,iJ,)(f) (d;q,n- 1) - %xt(.d;k,fi - l)<j) (d; q, fi) 

Here, we have used the common differentiation rules for Bessel and Hankel functions (see 
e.g. |54| , |3~TD . With the help of the Jost function we can now calculate the phase shift and the 
partial wave amplitude. We find 

J)(k ) - ' — - — - Q J n( kd )Ji*-i(q d ) ~ kJ tl _ 1 (kd)J fl (qd) 
^ 1 ,fl) ~ l F + F~ qY ll {kd)J^ 1 {qd)-kY^ l {kd)J tl {qd) 



f(k ) = 1. - F - ~ 2 g^M^zlM ~ kJ^i(kd)J^qd) 
,fi VmkF Vrik qH^\kd)J^ l {pd)-kH^\{kd)J t ,{qd) 

Inserting the last expression into the formula (|V.21| ) for the partial wave decomposition of the 
scattering amplitude we can calculate the differential cross sections for various parameters 
a. Note that for small momenta we can make use of the small energy behavior of the partial 
wave amplitude which is given by 

A;^0: f(k^) = 0(k 2 ^ 1 ). 
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Thus, when calculating the scattering amplitude in the partial wave decomposition the 
partial wave amplitudes for higher angular momentum hardly contribute to the infinite sum 
in ( | V . 2 1| ) . According to our discussion in section [TV] we have to add the Aharonov-Bohm 
scattering amplitude Jab m order to obtain the differential cross-section (see ( [IV. 8|) ). A 



numerical example is displayed in figure |j. For different values of the statistics parameter a 
one obtains very different angular distributions of the scattered particles. One observes the 
interpolation between the bosonic (a = 0) and the fermionic (a = 1) cross section at energy 
E = k 2 = 0.25. Note that the singularities at the scattering angles 6 = and 6 = tt for 
a ^ Z enter due to the long range nature of the statistics. In fact, according to ( |1V.£| ) they 



result from adding the Aharonov-Bohm amplitude given in ( [11.16| ). For fermions (a = 1) 



there is no scattering under an angle of 6 = tt/2 as is well known. The figure for a = 0.5 
shows that a similar effect can take place for fractional statistics. 

According to our discussion in section [V] the zeros of the Jost function F lie on the positive 
imaginary axis and determine the bound states. Thus, plotting the solutions of 

F(ik,fi)=0, k > 



we get so called (real) Regge trajectories (see e.g. f|3p displaying the /i dependence of the 
discrete spectrum of h^. Figure |5| provides a numerical example. In contrast to three dimen- 
sional scattering theory here every point on the Regge trajectories has a specific physical 
meaning because the statistics parameter a interpolates between the integer valued angular 
momentum channels. 



VIII. CONCLUSION 

We have demonstrated that it is possible to formulate nonrelativistic quantum scattering 
theories for two particles obeying anyon statistics. In particular we have proven the existence 
of the wave operators under certain assumptions concerning the interaction. For spherically 
symmetric potentials we gave criteria for completeness, i.e. unitarity of the S'-matrix. These 
criteria turned out to be the same as in three dimensional Schrodinger theory. It remains 
an open problem to show completeness in the general case in a way comparable to Enss' 
method. To generalize the latter, we believe it is crucial to obtain a better analytic control 
of the propagator. (There is an extensive literature on a similar problem, namely to obtain 
a closed form of the Aharonov-Bohm propagator see e.g. p"9|.) 

We extended the notion of differential cross-section to two dimensions and showed that 
in case of anyon statistics, the corresponding scattering amplitude consists of two parts. 
The first encodes the anyon statistics and resembles the Aharonov-Bohm amplitude, while 
the second is relevant to spectral analysis of the perturbed Hamiltonian. For spherically 
symmetric potentials we carried out this analysis by introducing Jost functions and showed 
how fractional statistics is equivalent to fractional angular momentum. We also showed that 
Levinson's theorem holds in the conventional case and gave the modifications necessary in 
the presence of a zero energy resonance. In the latter case we found that for positive angular 
momentum /i < 1 the statistics parameter - independently of the detailed form of the short 
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range potential - can be determined from the scattering phase, namely (compare ( |V.26| )) 



a 



— (5(k = 0, a) — 8(k = oo, a) j mod Z. 



We then applied our results to the following two examples. First we considered the 5— 
potential, where we extracted information regarding the point spectrum of the perturbed 
Hamiltonian from its resolvent and showed that the above relation between statistics param- 
eter and scattering phase is also valid in this context (compare (|VI.8| )). Second we considered 



the square well potential, where we obtained similar information from the corresponding Jost 
function. In both situations we numerically evaluated the differential cross-section for non- 
integer values of the statistics parameter. It then became evident that fractional statistics is 
fundamentally different from Bose and Fermi statistics: The statistical gauge forces produce 
singularities of the differential cross-section in forward direction, showing the long range 
nature of these forces. Due to this fact the angular dependence of the statistical interaction 
will invariably show up in the scattering data. In particular this might cause a differential 
cross-section which is constant for bosons or fermions, but which becomes angular dependent 
for intermediate statistics. Our discussion of the 5-potential provides an example. 

The results presented in this paper may also be relevant to the investigation of bulk proper- 
ties of anyon matter, via the relation between the virial coefficients and the scattering 
This we intend to discuss in a forthcoming publication. 
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APPENDIX A: THE BUNDLE THEORETIC FORMULATION 



In section [TJ we used the physical picture of particles carrying flux-tubes in order to 
describe the anyon model. Anyons, however, are but a special case in a more general class 
of particles, called plektons, obeying neither Bose nor Fermi statistics. In this appendix 
we give a short review of plektonic systems, i.e. systems of n particles whose statistics 
is determined by a finite dimensional unitary representation of the braid group B n . The 
appropriate mathematical tool to describe such particles is the concept of vector bundles. 
This will in particular provide the Hilbert space and a canonical "free" dynamics described 
by a certain Laplace operator. Our presentation will closely follow the one given in pjI)|.[TT 



where, however, emphasis was mainly given to relativistic formulations using momentum 
space considerations. For other references on the bundle theoretic formulation see e.g. 
We continue to use C to describe the one-particle configuration space. Let C xn denote the 
ra-fold product, viewed as the configuration space for n distinguishable particles. The set 
D n in C xn is the set of all points z = {z\, z 2 , . . . , z n ) with Zi = Zj for at least one pair of 
different indices. Any element a of the permutation g roup S n cicts in an obvious way on C xn 
via (za)i = z a ^. We define the configuration space of n identical particles in two dimensions 
to be 



"C := (C xn \ D n )/S n , 

with points in it denoted by z. Let ra C denote the universal covering space of ra C with points 
there being written as z and be 7r : n C — ► n C, z i— ► ir(z) = z the associated projection 
mapping. 

It is well-known that the fundamental group of n C, denoted by 7Ti( n C), is isomorphic 
to the braid group B n and that any element b G B n acts in a standard way from the right 
z i— > zb on the manifold n C. Furthermore, the universal covering n C i— > n C can be viewed 
as a principal fiber bundle over n C with structure group 7r 1 ( ra C) and fiber vr _1 (2;) for any 
z G n C Given any finite dimensional unitary representation b h- > p(b) of B n in a Hilbert 
space F with scalar product (• | •) there is an associated hermitian vector bundle. This 
vector bundle T with base space n C is given by n C X p ,B n F, which by definition is the set 
of orbits in n C x F under the following action of B n on this space 

b : (z, f) ^ (zb, p(r 1 )/) , z G "C, / G F (A.A1) 

T is a smooth fibered space with base n C and fibers isomorphic to F . On the fiber over 
z there is a natural scalar product denoted by (• | •) . Furthermore there is a canonical 
measure d/j(z) on n C inherited from the Lebesgue measure on C xn . This defines a scalar 
product on the space r c (^ r ) of smooth sections in T with compact support via 

<V> 10):= J &{z) | cj>{z)) z _dKi) ■ (A-A2) 

By £ 2 (jF) we denote the resulting Hilbert space completion. There is another Hilbert space 
equally well suited and canonically isomorphic to this space. Consider the set of maps 
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v& : n C i — > F which are smooth with 7r(supp ty) C ™C being compact and which satisfy the 
equivariance property 

V{zb) =p{b- 1 )V{z), Vze n C,beB n . (A.A3) 

For any two such functions $ and \1/ we therefore have 

(tf (z&) | 1>{zb)) = (tf (*) | $(5)) V z G n C, 6 G 5 n . (A.A4) 

Hence this expression depends onz = 7r(S) only and so the integral over the base space 
makes sense and we may define a scalar product by 

(tf | $) ;= y ($(5) | $(z))dfi(z). (A.A5) 

The resulting Hilbert space obtained again by completion is denoted by £^( n C, F) and there 
is a canonical isomorphism between £ 2 (JF) and £ 2 ("C, F) (see e.g. |4"0| , f4"Tp . Furthermore the 
canonical (flat) Levi-Civita connection on C xn induces a (flat) connection on n C which 
in turn defines a hermitian connection V on T. The associated Bochner or generalized 
Laplacean (see e.g. ||) A = Vo V on £ 2 (^-") is what defines a free Hamiltonian Hq = — A/2m 
for a system of n plektons, if m > is taken to be the mass of one particle. 



When the unitary representation p of B n is one-dimensional one speaks of anyons. Any 
such representation is obviously abelian and can be shown to be of the form bk i— > exp(ian) 
in terms of the generators bi, . . . 6 n _i G -B n for a fixed a G [0, 2). This follows easily from 
the observation that all the fr^'s are conjugate to each other. We denote the resulting line 
bundle by T a . Another consequence of F being one-dimensional is that all anyonic line 
bundles are actually trivial. This observation was first made by J. S. DowkerB based on 
V. I. Arnold's result that H 2 ( n C, Z) = oi and the classification theorem of Cartan, Kostant, 
Souriau and Isham. It was rediscovered by M. GaberdielH and is implicitly contained in [59 



and 17. Another proof is given in 40,41 



Let us now consider in more detail the Hilbert spaces constructed above when n — 2. We 
first introduce relative coordinates by considering the following transformation of C x2 : 

(zi, z 2 ) i — ► (^ 1 (z 1 + zi), z x - z 2 ) = (z cen , z re i) (A.A6) 

The transposition in S 2 — Z 2 obviously maps (z cen ,z re[ ) into (z cen , — z re {). Therefore 2 C is 
diffeomorphic to C x C*/Z2. The first factor is the configuration space for the center of mass 
motion. The corresponding quantum mechanical discussion is analogous to ordinary multi 
particle systems, since it is not affected by the statistics. Therefore we will concentrate on 
the second factor and the associated quantum mechanical description. In the same way as 
R is the universal covering space of S 1 , the universal covering space of C*/Z 2 = M + x S 1 is 
given by M + x K with the projection mapping 

7T : M+ x R — , with (r, 9) i — > (r, e 2id ) . (A.A7) 
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The action of B2 on the universal covering is given in terms of its only generator b\ as 

(r, 0) 1 — ► (r,0)&i := (r,9 + ir) . (A.A8) 

Choosing the representation h\ exp(— ma) for a fixed a £ [0, 1] we are prepared to 
construct the line bundle T a and the Hilbert space C 2 (T a ) or equivalently £L( n C,F). The 
fact that the bundle T a is trivial is now reflected in the possibility to "pull-down" the theory 
from either of those C 2 spaces onto the space of square integrable functions on n C itself, 
denoted by £ 2 (™C). Explicitly there is the following unitary map from £ 2 ( n C) onto C^ q ( n C) 
given as 

^£ 2 (T) — >^e£ 2 eq ( n C) with = e i7ra ' 0(5) • V (tt(2)) , 

where 0(5) is a real valued, continuous function of 5 induced by the polar angle. With the 
help of this unitary mapping, we can pull-down the Bochner Laplacean from C 2 (T a ) to 
obtain an expression for the free Hamiltonian on £ 2 ("C). The result of this procedure is 
precisely the Hamiltonian (O) we introduced in section [TT] on physical grounds. 

Let it be noted that a was introduced as a parameter fixing the representation B 2 3 
bi 1 — > exp(— irca) and is hence determined up to an additive even integer, only. Therefore 
we loose no generality if we restrict a to the interval [0,2). In fact, we have Ho(a + 2) = 
exp(— i26)Ho(a) exp(i26), which is a gauge transformation reflecting the arbitrariness we 
have. On the other hand, time inversion causes a to change sign, i.e. if T denotes the 
anti-unitary operator of time inversion, we have H (—a) = TH (a)T. However H (— a) 
is in turn equivalent to H (2 — a) and consequently the anyon models with a £ [0, 1] are 
connected to those with a £ [1,2] by time reversal. This justifies the restriction of a to the 
interval [0,1]. 

Furthermore we remark that £ 2 (C*/Z 2 , rdrdO) is unitary equivalent to the square in- 
tegrable symmetric functions on the punctured plane P + C 2 (C*, \rdrdff) and therefore the 
case a = corresponds to the bosonic case (since H (a = 0) = —A). If a = 1 our model is 
equivalent to H (a = 1) acting on P + £ 2 (C*), which is in turn unitary equivalent to —A on 
P_£ 2 (C*, ^rdrd9), the Hilbert space of square integrable, anti-symmetric functions on the 
punctured plane. Hence we call anyons obeying statistics corresponding to a = 1 fermions. 

There is also a geometric description of the universal covering space essentially as the Rie- 
mann surface of the function log(,z), which is useful for the understanding of the Hilbert 
space constructions made above and which goes as follows. Take an infinite set of upper and 
lower half planes in C and glue them together alternatively along the positive and negative 
axes respectively. In fact, consider an element / in C^° a (M \ {0}). Then it is easy to see 
that this / uniquely defines a function \1/ in C^° eq ( n C, F a ) and conversely any function \1/ in 
Cq^^C, F a ) defines a unique element / in C^(H\{0}). Furthermore this correspondence is 
linear and the scalar products are the same under this correspondence. This establishes the 
desired isomorphism between TC a and £ 2 g ( n C, F a ) and hence also C 2 (T a ). This isomorphism 
also extends to the various Laplace operators. Indeed, note first that locally n C looks like 
C xn , so there is a natural Laplace operator A on C°°( n C), the space of smooth functions on 
n C. This operator is easily seen to map ( n C, F a ) into itself and hence defines an opera- 
tor on C 2 ( n C,F a ). Furthermore this operator also corresponds to the Laplace operator on 
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Co°(IHI\{0}) under the above correspondence / «-> iff. Finally it is easy to see that under the 
canonical isomorphism between Cl q ( n C, F a ) and £}[Tq) referred to above A corresponds to 
the Bochner Laplacean A = V o V. 



APPENDIX B: SOME ESTIMATES 

This appendix is devoted to a proof of the lemmas [III.2| and [III. 3 . 



Proof of lemma [7?Q ; We start with a proof of lemma |111.2| . For the moment let < e < 1 
be arbitrary. Since \e lx — 1| < min(2, | x |), we have \e lx — 1| < 2 \x\ e for all i6l. This 
gives 



\Up,x)-Uo,x)\ < 



dy ( e ipcoshy -1) - 



-ay 



< 2p e / dy cosh 6 y 
Jo 



-ay 



-2y-2i x 



+ 



?ay 



I — e -2y-2i x I _ e 2y-2i x 

To proceed further, we split the integral into a part from to 1 and the rest. This gives 

\I a (p, X )-Ia(0,x)\ <2p e (cosh(l)-/ 1 ( X ) +J 2 ( X )) . 
For Ii we have an expression involving the sum of three integrals Ix % and Jj 3, given by 

+ 





e~ ay - 1 




e ay - 1 


/ 

to 


+ [ dy 
Jo 


1 _ e -2y-2i X 


1 _ e 2y-2i x 



1 

I dy 
'0 



1 



1 _ e -2y-2i X I _ e 2y-2i X 



For I2 we have the expression 

POO 

h{x) = dy cosh e y 



-ay 



I — e -2y-2i X 



+ 



1 _ e 2y-2i x 



The aim is to estimate these quantities for x £ [ — vr, +7r]. To estimate Ii,i(x) we use the 
following little estimates: 



|1_ e -*y\ < 2y 

| 1 _ e -2y-2i X > h_ e -^| 



for y > and a e (0, 1] 
for y eR: 



(B.B1) 
(B.B2) 



This gives 



h,i(x)< J o dy- 



ay 



-2.y 



< c . 



h,2(x) is estimated similarly, if one replaces the estimate ( |B.B1| ) by 

< e ay - 1 < 2e 2 y , for < y < 1 and a e (0, 1) 
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To estimate ii,3(x) f° r X £ l~ 7V i + 71 "] we first observe that for min( \x\ , 7r 
it is obviously bounded. Hence it suffices to estimate the three remaining cases 

3 



X^ + X) > f 
< — 

4' 



f 71 " < X < 7T or — 7T < x < — j7r. We only consider the first case, since the other two 



\X\ 



< ~7T and add and subtract 



cases may be discussed with similar arguments. So let 

(±2y + 2ix)~ 1 in the integrand of /^(x). This way we can estimate I\ t 3 by a sum of three 
integrals, which we will denote by Jj. 



Ii, 3 (x) <h{x)+Ux)+Ux) 



dy 



1 



1 



1 - e^v-^x 2y + 2ix 



dy 



1 



+ 



dy 



1 _ e 2y-2i X _2y + 2«X 
f 1 



2y + 2%x -2y + 2i X 



To estimate Ji(x) an d ^(x) we note that the function 

1 1 

C 



is obviously analytic in G C | |Im(^)| < |} except possibly at the origin ( = 0. However, 
G(C = 0) = 1 and hence by Riemann's theorem, G(() is an analytic function in that domain. 
In particular G(() is bounded on every compact subset and this gives the boundedness of 
h{x) and h(x) for x G [— -r^, j 7r]. hix) can De estimated as follows 



Ux) < / dy 



\X\ 



x +y 



< 



dt] 



1 + V 2 



< oo 



This concludes the estimate for h(x) arid it remains to estimate ^(x)- By ( |B.B2| ) we have 



Ux) < dy2e 



-ay 



1 - e~ 2 y 

Now this integral is finite whenever e < min(a, 1 — a). 



+ 



,2y 



□ 



Proof of lemma \uT. j| : We will consider the third partial derivatives w.r.t. E only, since this 
is actually the case we will need and since the other cases may be discussed similarly. By 
interchanging integration and differentiation we formally have 



d" 



3E r 



\E,E>) 



d" 



dE r 



4>a;m,E{ Z ) V ( Z ) <t> a;m' ,E' {%) djj(z 



(B.B3) 



for n = 1,2,3. This is permitted provided the integrand is a measurable function in z, E 
and E', bounded in the sense of the modulus by a function which is in C 1 w.r.t z and with 
uniform bounds in E and E' in compact sets in (0, oo). To prove this claim it obviously 
suffices to replace ^ by J| with E = k 2 in ( p.B3| ). By the explicit form of (f) a -m,E given 
by (|IlD, for fixed a, m and m' we therefore have to estimate the product 



\d n k J,(kr)\ ■ \ J,(k'r)\ ■ |y(r,, 



for n = 1, 2, 3 
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Using the formula to be found i 

^ J °(kr) = \ { J °-ii kr ) ~ Ja+i(kr) ) (B.B4) 

it follows by iteration that it suffices to estimate 

r n \J, +l (kr)\ ■ \J^k'r)\ ■ \V(r,9)\ 2 (B.B5) 

for —n < I < n with < n < 3. Let the compact set in question be given by Ei < E , E' < E 2 
such that ki < k, k! < k 2 . Using the well known estimates 

a) \Ja(kr)\ < C(a) for \kr\ < 1 

b) \J<r{kr)\ <C(a)(kr)~^ for \kr\ > 1 

we see that (|B.B5|) with k\ < k, k' < k 2 is bounded by 

C(ji, //) • r" \V(r,6) | 2 < C(ji, //') • max(l, A;^ 1 ) 3 \V(r, 9) | 2 

for r <k 2 l and by 

C(ji, //) ■ T n ~ x k~ x x \V(r, 9) | 2 < C(ji, //) • max(fc 2 , k?) • r 2 |U(r, 0) | 2 

for r > &7 . Obviously ( B.B5| ) is bounded for r in the interval (k^ 1 , k^ 1 ) uniformly for 
k\ < fc, A;' < ^2- This proves the claim and by our previous remark this concludes the 
proof. □ 

APPENDIX C: BOUNDS FOR JOST FUNCTIONS 



In this section we will prove the estimates ( |V . 1 1 b| ) , ( | V . 1 7| ) , ( |V.lla| ) and (|V.16| ) for the 



functions <fi and an d their power series coefficients <p n and x ±,n - Also we will show that 
the analyticity properties for <p and x ± a l so extend to their derivatives with respect to r. 
We start with two preparations. Observe first that (|V.llrj| ) and ( |V.17|) hold when V = 0, 



i.e. for 6 = <fi and = Xo respectively. Indeed for <fio ( |V ■ 1 1 b| ) follows from the relations 
(see e.g.S'0) 

cosiz - f (2tt+ 1)) 

J^z) > C(//)^ and J^z) 



z^O \ z \-+oo ^ T^Z 



For Xo the equation ( |V.17| ) follows from the following asymptotic behavior 
H±(z) — ► C{p)# and H±{z) 



z^O ' z— >oo + / —Z 

V 2 



Secondly for g (see (|V.7| )) one establishes analogously 

0<r'<r : \g (r,r'; M)l < C(ji) e IM*)IC~ ') (i^)^ , (G.C1) 

< r < r' : |<7 (r,r'; < CQi) eU-WI-M^) j^)^ f^)^" . (C.C2) 
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Now (|V.lla|) is proved by induction on n as follows. Define A n by 



l+|fc|r 



\<MrM\ = e |Im(fe)|r (t^ 

and it suffices to show that 



A n (r;k,fj,)<\k\^ 



nl 



dr 



A n (r; k,fi) 



, r'\V(r')\ 
1 + \k\r' 



where C(fi) is the maximum of the C([a) in ([V.llbD and the one in flC.Cl|) . By the prepara- 
tory remarks (|V.lla|) holds for n = 0, i.e. for 0o = <p n =o- To perform the induction step, by 
construction we have 



<f>n(r;k,fj) = / dr'g (r,r'; fc,/x)V(r')0„_i(r-'; k,/i), 
Jo 

which gives the inequality 

Mr; k, fi) < C{im) f dr' ^£2} A n _ x {r'; k, fi), 
Jo 1 + l ft l r 

from which the induction step follows easily. This completes the proof of ( |V.lla|) and 
( |V . 1 1 b| ) follows from it by summation. 

The proof of the bounds (|V.16|) and (|V.17|) is similar and we will consider the case x~ only. 
Now write 



such that it suffices to prove 



-/X+7 



B n (r; k,/i), 



B n (r; k,fi) < 



C{fi) n+1 



TV. 



\k\ 



°V *W)I N " 
1 + \k\r' 



By the preparatory remarks, this inequality is valid for n — 0. To make the induction step, 
we note that by construction 



Hence one has the estimate 



dr'g (r, r'; k, //) X n -\{ r '\ k , I- 1 )- 



F 00 r'\V(r')\ 
B n (r;k,fi)<C(n) / dr . B w _i(r ; fc, //), 

J r 1 + \k\r 

from which the induction step easily follows. This completes the proof of (|V.16|) and (|V.17|) 
follows by summation. 
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APPENDIX D: THE RESOLVENT OF THE 5 POTENTIAL 



Proof of theorem \ VI. For any rj G £ 2 (M + ) define 



V 



using the integral kernel given by ( |V1.2| ). It is easily seen that is well-defined. In order 
to prove the theorem we shall show that (|VI.2|) formally defines a resolvent, i.e. we have 



d 2 " 2 -±> 
■- — V -k ip = rj 



(D.DI) 



That this "resolvent property" is formally satisfied, can easily be verified with the aid of the 
well-known formulae for the first derivatives of Bessel functions. In particular one uses the 
relations: 

C a = — C a — C a +i Cq+i = C a - — Cq,+i for C a = J a or H a . 



z 



z 



Furthermore we shall show that ip satisfies the boundary condition in ( [VI. lj ), i.e. 



limW(V>«, a ,V)(0 = . 

rj.0 



(D.D2) 



For convenience we introduce the following notation: 



J>(r) := 



I<(r) :— j dr' Vr' J a (kr')rj(r') 
dr' V? H^\kr')rj{r' 
ip can now be cast into the form 
ip{r) = \fr 



POO 

I ■= dr' Wr'H^\kr')r]{r') 
Jo 



yffW(Ar)/<(r) + yJ Q (fcr)l>(r) - A(k, a; s)H^\kr)I 



(D.D3) 



To verify the boundary condition ( p.D2|) we note that 



VVs( r ) = — " r2 a + s-r 2+a with 
2a 



s = s 



2a 



giving 



- + sa ) r ~ 



and hence W{tp a , s , ip){r) = B\r) ■ r 1_Q + 5r 1+c ^j + B(r) Q ■ r" Q - asr°^j 



where B(r) denotes the quantity in the square brackets in (|D.D3| ). In order to take the 
limit, we remark that z~ a J a (z) is an analytic function for any value of a 6 1. From its 
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power series expansion at z = (see e.g£3) one obtains the following asymptotic relations 
when a > 0: 



limr a J a (kr) 
lim r a Y a (kr) 

r|0 



lim r a Y_ a (kr) 

r|0 



L a (k) := 
la{k) :-- 

l-a{k) 



1 



2y r(l + a) 
2\ Q 1 1 



J sin 7TQ; T(l — a) 
k\~ a 1 1 



2 / tan7ra T(l — a) 



For the definition of the integral kernel ( VI. 2 ) we restricted the Bessel functions to the sheet 
given by < arg(A;) < tt. Since the irrational powers of k in the above expressions stem 
from the asymptotic behavior of the Bessel functions, they have to be evaluated on the same 
sheet. Bearing this in mind, we are prepared to take the limit r — > + , which will make 
most of the terms in ( p.D3| ) disappear. Since |/<(r)| ~ C ■ r a+1 and i> ~ / for r < 1, the 
surviving terms are: 



lim W(ip ayS , ip) (r) = I ■ lim 

r|0 rjO 



ITY 

Y 



r~ a J a {kr) + s-k- A{k, a; s)r 1+a H^ +l {kr) 



Now, by a standard theorem for Bessel functions, we have H^}.i(kr) — ^Ha\kr) 
—Hal\{kr). The boundary condition (p.D2j) therefore implies 
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A(k,a; s) = —i a {k) 



i a (k) - 2isa ■ 7„ + — 7q-i 
2a 



-i 



Substituting s = s — and the expressions for L a , 7 a and 7 a _i we see that A(k, a; s) is 



given by ( |V1.3| ), thus completing the proof of theorem |V1.2 



□ 
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FIG. 1. The diagram shows the location of the square root of the bound state 
energy kb on the imaginary axis as a function of the coupling parameter s for ten 
different, equispaced values of the statistical parameter a. "bosons" labels the graph 
corresponding to a = 0, "semions" corresponds to a = |, the graph intersecting the 
s-axis at 1 and "fermions" corresponds to a = 1. The latter case does not correspond 
to a bound state. 



46 




FIG. 2. The above figures show the "normalized" differential cross-section -r-y — % 

° sm z na at) 

as a function of the scattering angle 9 in units of it. The normalization factor is chosen 
such that the asymptote in the limit s — > — oo (i.e. when the interaction is turned 
off, giving the Friedrich's extension of ho jCe ) becomes independent of a. The latter is 
shown in the first plot (top left) corresponds to Aharonov-Bohm scattering. The other 
five plots show a family of graphs corresponding to the following coupling strengths 
s = — 10, —8, ... , 0, k = 1 and the given values for a. For small values of the graphs 
for different s are ordered from left to right starting with s = 0. 
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FIG. 3. Similar to figure y the above plots show the differential cross-section 
sin^Tra W as a f unc ti° n of the scattering angle 6 in units of ir. But now only posi- 
tive s values are considered. The first plot (top left) shows the limiting case s — > oo 
corresponding again to Aharonov-Bohm scattering. The other five plots display a fam- 
ily of graphs with coupling strengths s = 0, 2, . . . , 10, k = 1 and the given values for a. 
Here the ordering of graphs is different from figure ||. For small values of 6 the leftmost 
graph always corresponds to s = but the ordering is reversed once at a certain value 
s = s(a). 
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FIG. 4. Displayed are the differential cross sections per unit length (y-axis) in 
dependence of the scattering angle 9 in units of ir (x-axis). The first figure shows the 
symmetrized Aharonov-Bohm cross-section divided by sin 2 ira for non-integral values of 
a. The other graphs depict the differential cross-section when the statistical interaction 
as well as the square well potential are present. The depth of the square well is chosen 
to be Vq = 25 and the radius is set to d = 1. 
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FIG. 5. The graphics shows the zeroes of the Jost function F in dependence of 
the momentum k = in, k > and the angular momentum \x for the square well po- 
tential. The line broadening is due to numerical error. The chosen numerical values 
are d = 1, Vo = 25. One observes two bound states the first of which disappears at 
\x = 0.674 and the second at \x = 2.893. 
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